Bose-Einstein condensation under external 

conditions 



Klaus Kirsten [flj 
University of Leipzig 
Institute of Theoretical Physics 
Augustusplatz 10, 04109 Leipzig 

and 

David J. Toms |f 
Department of Physics 
University of Newcastle Upon Tyne, 
Newcastle Upon Tyne, U. K. NE1 7RU 

February 1, 2008 

Abstract 

We discuss the phenomenon of Bose-Einstein condensation under 
general external conditions using connections between partition sums 
and the heat-equation. Thermodynamical quantities like the critical 
temperature are given in terms of the heat-kernel coefficients of the 
associated Schrodinger equation. The general approach is applied to 
situations where the gas is confined by arbitrary potentials or by boxes 
of arbitrary shape. 
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One of the most striking developments of the last few years has been the 
experimental observation of Bose-Einstein condensation in very cold gases of 
rubidium jE), lithium and sodium ||. The systems are very dilute and 
as a first approximation would be expected to be well described by a boson 
gas model with no interactions among the atoms. The atoms are confined 
in complicated magnetic traps which can be modeled by harmonic oscillator 
potentials. Nowadays one can say that the thermodynamic properties of the 
confined gas are very well understood, including the effects that the finite 
number of particles has on the system ||, [5| ||, [F], |], |] as well as the influence 
of a weak interaction among the atoms [|nj. The above mentioned results 
have all been obtained in the grand canonical approach, which, although 
inadequate for the calculation of, for example, flucuations of particle numbers 



[IT], [12], [LJI |14]] is still reliable for thermodynamical properties such as the 
critical temperature or the ground state occupation number. 

Although recently most interest is in Bose-Einstein condensation in mag- 



netic traps also different confining potentials [15| or gases confined inside a 



rectangular box |16[ have been considered and specific differences have been 
stated. A formal discussion of boundary effectis has been given in [|17|]. All 
this research may be called analysis of Bose-Einstein condensation under 
external conditions. 

In this article we want to provide an approach which enables one to deal 
with any kind of external conditions (including the most relevant anisotropic 
harmonic oscillator potential). The methods we are going to employ are the 
so-called heat-kernel techniques extensively used in finite temperature rela- 
tivistic quantum field theory starting with the work of Dowker and Kennedy 



|T8| . (See also PD|| .) in non-relativistic theories these techniques are however 
nearly unemployed and we want to show in the present article that they can 
also be used here in a very effective way. 

Let us consider a system of N non-interacting bosons under external 
conditions such that the single particle state energies En are determined by 
the Schrodinger equation 

h 2 

- —A4> N (x) + V(x)(j) N (x) = E N (p N (x) , (1) 
2m 

where V(x) is the confining potential. For example in [15| the confining 
potential 

V{x) = jZef\x i r (2) 
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has been considered in the thermodynamic limit. If in addition the gas is 
confined by a box Ai one has to impose suitable boundary conditions at the 
boundary dM. of the box. Without specifying the spectrum En (as a rule it 
will not be known explicitly) the grand canonical potential may be given in 
the form 



OO 1 

9 = <?o + £ £-e-^"-"> . (3) 

n=l N U 

Here 

go = -d ln(l-ze-P Eo ) , 

is the contribution of the ground state with energy E , do is its degeneracy 
and z = exp(/9/x) is the fugacity. The prime in (|3|) indicates that the ground 
state contribution is to be omitted from the sum. 

For the evaluation of this kind of expression it is very effective to make 
use of the Mellin-Barnes integral representation pOf 

c+ioo 

e~ v = — [ da r{a)v~ a , (4) 
2m J 



C—lOO 



valid for dtv > and c 6 R , c> 0. This is easily proven by closing the 
contour to the left obtaining the power series expansion of exp(— v). As an 
immediate consequence of eqs. (|3]) and (|4]) the grand canonical partition sum 
may be casted into the form 



c+ioo 



q = % + ^n / ^rWr^'n-a^^jCla), (5) 

c—ioo 

with the polylogarithm 

oo x l 

^nW=E^ ( 6 ) 
1=1 1 

and the spectral zeta function 

aa)=Y,'(E N -E y a . (7) 

N 
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Here we introduced the variable fi c which denotes the value of the chemical 
potential at the transition temperature. For ideal gases this value is identical 
with the ground state energy, \i c = E . The parameter c in eq. fl5|) is given 
in a way such that all poles of the integrand lie to the left of the contour. 
Closing the contour to the left one picks up the right-most residues of ((a). 
Concerning their determination there are deep connections between the zeta 
function of an operator and its heat-kernel defined as 



K{t) = £' e -^- £ °). (8) 

The small-t behaviour 



N 



K(t)~J2ht- jk , (9) 

k=0 

is of particular relevance here, b^ are the so-called heat-kernel coefficients. 
We assume jk > jk+i, and the connection between the residues of ((a), which 
occur at a = jk, and the heat-kernel coefficients reads 

Res C(a = jk) = p^y. (10) 

Taking into consideration only the two right most poles of ((a) one arrives 
at 

q = q + p- j0 Li l+j0 (e-P^-ri) b + r jl Lh +jl (e^^) b 1 + ... . (11) 

For the determination of the critical temperature and the thermodynamical 
properties of the gas up to T c the above approximation can be shown to be 
a very good one. 

We thus have a representation of the grand canonical partition sum in 
terms of the heat-kernel coefficients of the relevant Schrodinger operator. 
What have we gained with this representation? This becomes clear by notic- 
ing that the heat-kernel coeficients bo, b\ are known or easily determined for 
a very large class of Schrodinger operators including the situation where the 
gas is confined by arbitrary potentials or boxes of arbitrary shape. This in- 
cludes physical situations which are of relevance in recent experiments and 
also those which might be of relevance in future ones. 



4 



Let us exemplify the power of the approach by considering further the 
particle number and the critical temperature. The particle number is 



N 




\T,V 



= N + (3- jo Li jo (e-^-")) 6 + P~ h Li h (eT^-A) b x + ... (12) 

In order to avoid distinguishing all the different cases, let us assume from 
now on that j\ > 1. In this case, the critical temperature 1/(3 C is determined 
by 



and reads 



T c = T 1 - 



Cj?C/'i)6i 



(13) 



(14) 



with 



T n = - 



1 / N 



(15) 



the critical temperature in the bulk limit. Let us stress that eq. ( 14]) con- 
tains the influence that the finite number iV of particles have, the details 
being encoded in the exponents ji and the coefficients 6j. Proceeding with 
the energy and the specific heat, all thermodynamical quantities are easily 
displayed in terms of the heat-kernel coefficients. If ji < 1 different results 



are found due to the behaviour of the polylogarithm [21 



The presented technique allows us to calculate thermodynamical proper- 
ties by treating the sums over the discrete energy levels. Another possible 
way to do this analysis is to approximate the sums by integrals. A crucial 
feature in obtaining a reliable approximation is to use an appropriate density 
of states 0. By using a refinement of Karamatas' theorem one can show 
that the use of the density 



p{E) 



bo 



61 



-E jl 



r(jo) m) 

is completely equivalent to the analysis presented above. 
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A slightly different approach results by the use of the effective fugacity 
(for the anisotropic harmonic oscillator see 0) 

z eff = ze-™ 

where E\ is the first excited level with, lets say, degeneracy d%. Whereas in 
the previous calculation only the groundstate has been treated separately, 
we now separate the groundstate and the first excited level to find 

q = q + dxLi^Zeff) + (3~ ]0 Li 1+jo (z e}f )b + f3~ n Li 1+jl (z eff )bi + ... 

This expansion is a very good approximation even below the critical tem- 
perature (for the harmonic oscillator see J7|). Also here thermodynamical 
properties are obtained in the same manner as before. 

Let us now use our general results to consider the special situations men- 
tioned already several times. 

Example 1: Gas inside a <i-dimensional box of arbitrary shape. 

In the notation of equation ([!]) the situation is described by V(x) =0 and 
4>N( x )\xedM = 0- (Other boundary conditions can be treated as well.) For 
this example the relevant heat-kernel coefficients are known (see f.e. ||23|| ) 
and read 

b = (4ir)- d/2 (2m/h 2 ) d/2 vo\(M) 

bj_ = -(4Tty d / 2 (2m/h 2 )( d - 1)/2 (V^/2)vo\(dM) 

together with j = d/2, j\ = (d — l)/2. Here vo\(M.) and vo\(dA4) denote 
the finite volume of the box and its boundary respectively. For the critical 
temperature this means 

Tc = Tq f 1+ 1 vol{dM) Ca((d-l)/2) 1 ] 



2d (vol(M)Y d - i y d C R (d/2)( d - 1 )/ d N l / d 
where 

h 2 ( N \ 2/d 



T 



2nmk \vol(M)( R (d/2) 



) 



The rectangular three dimensional box has been analysed in detail in [TJ 
using a density of states approach. 



6 



Example 2: c/-dimensional isotropic harmonic oscillator as confining 
potential. 

Let us now choose the confining potential 



V(x) = -Tiujm^2, 



X, 



In this case the coefficients are found by just calculating with geometric series 
0. The result is b = {fiui)~ d , b\ = (hdu/2)b and j — d, ji — d — 1 and so 

t-tU- C^-i) l \ 

c °\ 2Ca{dy d - l )l d NV d \ 



with 



N \ 1/d 



(see also f.e. [0). 

Example 3: 3-dimensional anisotropic harmonic osciallator as confining 
potential. 

In the anisotropic case the relevant potential reads 
V(x) = -%m{u)ix\ + lo 2 x\ + oj%x%). 

With the same comments as above one gets here jo — 3, ji = 2, and 
1 1 ( 1 1 1 

b o = 73 ; b i = 7^2 [ 1 1 

n Uiuj 2 uj 3 2n \U\U2 ujiuj^ uj 2 ^z 

For the critical temperature this means 

T =rjl ryN- 1 / 3 ] 

c °l 3Cfl(3)"/3 77V J 



with 



and 



7 = ^(ujiu 2 uj 3 ) 2/i 



1 1 1 
— + + — 



UJ1UJ2 CU1CU3 UJ 2 UJ 3 



T = ft(a>ia> 2 a> 3 ) 1/3 
7 



N \ 1/3 



These results have also been obtained using an approach based on the Euler- 
Maclaurin formula ]7) and based on a density of states approach [[J. 

Analogous formulas for the <i-dimensional anisotropic harmonic oscialla- 
tor can be found in an easy manner. 

Example 4: Arbitrary power law confining potential in 3 dimensions. 

Let us finally consider the potential 



v 
1 1 •> 



where in order to simplify comparison with the harmonic oscillator calcula- 
tion we use the notation 

m, 2-P £+2 

—ll ~<" Jj -<• . 



Using the results of [25] for the resumed heat-kernel, one gets for this case 

_ sr 3 (i/ P ) i 

° vr 3 /V (fto;)f(l+2/l0 

and b\ = E bo, together with j = | + - and ji = jo — 1. The resultm; 
critical temperature for this example is then easily found to be 



(E /MCfl(l/2 + 3/p) / 2r(l/p) \ 



2p_ 
P+2 I 



3(p+2) 



T r . — T < 1 p+6 I 7= I 7T 

(3/2 + 3/j9)Ci?(3/2 + 3/p)WT V V*P J N — 
The bulk critical temperature here is 



2p 2p 



JV \ 3(p+2) 



y/np y 



T = Tio; 



(for further bulk properties see |I5|). 

In summary we have shown, that our approach enables us to deal with 
complicated external conditions in an efficient way. Of course we do not 
claim that the results presented in the examples cannot be obtained by other 
means; in fact, most of the results for T c have been obtained differently. How- 
ever, the unified treatment of many different situations seems very elegant to 
us. Having at hand this general approach given an ideal Bose gas under ar- 
bitrary external conditions one only has to determine the leading heat-kernel 
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coefficients of the associated Schrodinger operator (or even only take them 
from the literature) in order to know the whole thermodynamical behaviour 
up to or below the critical temperature. Without technical complications our 
approach includes the effect of a finite number of particles. 

These ideas are also applicable to the canonical ensemble. Furthermore, 
due to the connections between the microcanonical and the grand canonical 



ensemble [13[] it seems possible to develop this approach also for the micro- 



canonical ensemble. 
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